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Abstract 

Properties of confined mesoscopic systems have been extensively studied nu- 
merically over recent years. We discuss an analytical approach to the study 
of finite rotating fermionic systems in two dimension. We first construct the 
energy functional for a finite fermionic system within the Thomas-Fermi ap- 
proximation in two dimensions. We show that for specific interactions the 
problem may be exactly solved. We derive analytical expressions for the den- 
sity, the critical size as well as the ground state energy of such systems in a 
given angular momentum sector. 
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I. INTRODUCTION 



The observation of Bose- Einstein Condensation(BEC) in extremely cold dilute gas of 
atoms has generated a lot of activity both on experimental as well as theoretical front on the 
properties of trapped dilute gases. There is now a renewed focus on the properties of trapped 
dilute gas of fermionic atoms at low temperatures. A recent effort in this direction is the 
experimental observation of quantum degeneracy in a dilute gas of trapped fermionic atoms 
||T|]. Several theoretical papers have studies the properties of trapped dilute gas of fermionic 
atoms. Butts and Rokshar |^ have studied the momentum and spatial distribution of the 
noninteracting system in the Thomas- Fermi approximation. The ground state properties and 
the addition energy spectra of a two-dimensional interacting fermi system has been studied 
by Sinha et al 0. Recently collective excitations of the system in the normal phase [§] and 
in the superfiuid phase have also been investigated. The thermodynamic properties of 
fermi gases including rotations has been analysed by Salasnich et al 

Especially relevant to this paper is the recent interest in the properties of rotating bose 
condensates and the experimental observation of vortices in stirred BEC of ^^Rb . The 
properties of ground and low excited states of a rotating, weakly interacting Bose-Einstein 
condensate in a harmonic trap has been investigated recently in various limits [^|10|- Ho 
and Ciobanu ||T^ have recently examined theoretically the nature the ground state of a 
rotating trapped fermi gas in two and three dimensions in the noninteracting limit. They 
show that the density profile acquires features reflecting the underlying Landau level like 
energy spectrum. Properties of the spectrum of such a system had been earlier investigated 
by Bhaduri et al [Q. 

In this paper we investigate a system of rotating fermions in two dimensions confined 
within a parabolic potential with some model interactions. We show that within the Thomas- 
Fermi approximation the system may be solved exactly and discuss the consequent results. 



We follow closely the formalism developed by Gallego and Das Gupta [|T3[ who developed 
the Thomas-Fermi method for rotating nuclei. The attractive feature of Thomas-Fermi 
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method is the ease with which non-trivial many-body solutions can be obtained even when 
they are interacting. As we show later, in some cases the results may be obtained purely 
by analytical methods. While this is not a full list of the exactly solvable models, they 
are chosen for their relevance to physical systems. The paper is organised as follows. In 
section II we outline the derivation of the Thomas Fermi energy for a system of interacting 
rotating fermions in parabolic confinement. In section III we discuss the results for three 
specific type of interactions between the fermions. We obtain analytic expression and present 
numerical results for the interaction energy, the rotational energy and the spatial density in 
a given angular momentum sector for a given number of particles. The last section contains 
a discussion of the results. 



II. THOMAS FERMI ENERGY FUNCTIONAL 

We derive the Tho mas- Fermi (TF) energy functional for a confined two dimensional ro- 
tating fermi gas starting from the microscopic Hamiltonian. We shall follow the approach 



outlined by Gallego and Das gupta [|T3[ who derived the TF density functional for rotating 
nuclei, but applied to a two dimensional system. 

The microscopic Hamiltonian of a two-dimensional rotating fermi gas confined in a har- 
monic potential is given by. 



N „2 



2„2 



0' i 



1=1 i<j 

- ^mu;[ri X Pi], (1) 

where V{ri,rj) denotes the two body interaction between fermions, ujq and uj denote the 
confinement and rotational frequencies respectively. 

We want to minimise the energy of this system subject to the constraints 



drp{r) = N, (2) 
j dfdpf{f,p)(fx p) = L, (3) 
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where N and L are total particle number and the total Angular momentum, f{f,p) is the 
semiclassical phase space density which we shall define later. The configuration space density 
p is obtained by integrating the phase space density over the momentum space. 

Using Lagrange multiplier /i and cu for each constraint respectively and regrouping the 
momentum dependent term we obtain the following expression for the energy functional 

dfdpf{r, p) [— - uj{xpy - yp^)] 
+ j df\y{r) - iJLp{r)] + iiN + ujL (4) 

where the mean-field one body interaction V{r) is 

V{T)^hmujy + VH{,r). (5) 

The Hartree term is given by, 

VH{f) = J d'^rV{f,f')p{f'). (6) 

We ignore the exchange correction as it is in general subdominant compared to the Hartree 
term. 

The semiclassical phase space density is given by 
where we have accounted for a factor of 2 due to spin degeneracy and 

denotes the energy density in phase space. We may rewrite this by completing the square 

as, 

e(p, r) = — [{px + muoy)'^ + {py — muux)'^] + V{r) muj'^r'^. (9) 

The shifted momentum simply indicates that the center of the Fermi sphere at any given 
point f in the rotating frame is displaced from the usual p — 0. This shift is of no relevance 



within the classical TF approximation except for the appearance of the extra centrifugal 
term in the energy functional. 

The TF energy functional is therefore given by, 



E[p\ = f d^r [^—^ + -mnr^ + Vh- /i]p(f) +fiN + ujL (10) 
J 2m 2 



where 



fi2 = cul-iu^ = Ctu^ (11) 

denotes the effective frequency in the presence of rotations. The TF equation for the spatial 
density is obtained by a variation of the energy functional, namely, 

lll^ + ^rnn^r' + Vh = (12) 
m 2 

To obtain the ground state density in TF approximation one has to solve the above equation 
self-consistently with the boundary condition that the density vanishes beyound the classical 
turning point Tq, that is 

p(r) = 0, r > tq. (13) 

Note that in the absence of rotations Q = ujq and the solution of the above equation describes 
the ground state in the absence of rotations which has been analysed in detail before 0. 
We may remark that in general the above equation is difficult to solve analytically except 
in specific cases which we shall discuss later. 

The angular momentum carried by the system is given by, 

L--ISr^,nrM...-y..)--^I^Mnr^ (14) 

which is simply the classical expression for the angular momentum. We may therefore 
identify the Lagranges multiplier u with the angular frequency. 

We may also define the related energy functionals. The free energy of the system is given 

by. 



F = E - fiN -uL (15) 

The TF equation (|T2|) may also be obtained just by varying the free energy. We also define 
the energy in the rotating frame through 

E = E' + uL. (16) 



where E is the total energy in the static frame. With the help of eq. ([T^) it can be shown 
that irrespective of the nature of the interaction E' is given by, 

E' = ^fiN + imfi^ J <fr r^p. (17) 

This concludes the basic TF formalism applied to rotating fermionic systems. In the fol- 
lowing section we solve the TF equation for specific cases. Note that it is necessary that we 
choose u < ujq as otherwise the system becomes unstable. 



III. SOME EXACTLY SOLVABLE MODELS 

We now discuss special cases of interacting fermi gases which may be solved exactly. 
These models are chosen not only to facilitate analytical calculations, but also because in 
some limits these approximate the effective electron-electron interaction in a two-dimensional 
fermionic system. We note that the effective interaction may be dependent on the device 
characteristics as well as the number of particles. As such we find that it is useful to work 
with different approximate forms of interactions. Throughout we assume the system is spin 
unpolarized. The partially and fully polarized system in the presence of magnetic field will 
be discussed in a subsequent paper. 



A. Contact Interaction 



The first and the simplest is the contact interaction. In general any short-range interac- 



tion may be written in terms of a gradient expansion [|1^]. The leading term in the expansion 
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is always the contact interaction of the form 

V{r) = K,5(r), 

where f denotes the relative coordinate between any two fermions. In a gradient expansion 
the factor Vq is related to the leading moment of the potential |T3|,|TB . 
The free-energy is then given by, 

F[p\ = J dh{^^ + ^mfirV(r) + yp'] - ^l j d^rp{r) (18) 

The TF equation for the spatial density then becomes 

+ ^]V^^, (19) 

m 2 

where g = which is dimensionless. Note that eq.(0) is identical in form to the case 
when there is no rotation. Only the effective frequency Q is different now. The solutions 
are therefore obtained easily by merely replacing the confinement frequency by the effective 
frequency. Since the density vanishes beyond a turning point tq, the chemical potential is 
related to the turning point by, 

= ^mQrl. (20) 

The turning point and hence the chemical potential are therefore determined by the total 
number of particles N which is fixed. 

Substituting it back in the TF equation gives the expression for density 

The total number of particles is then obtained by integrating the expression for density 
in eq.(pTl) 

N = rd\p[r) = (22) 
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and the total angular momentum is given by 

a rmujr pir) = tt- (23) 

In a system where both and L are fixed, the turning point tq and uj are determined by 
entirely by these two constraints. 

Using the expression for density and the total angular momentum L in eq.(|23|) we can 
easily evaluate the expression of the energy. Some straightforward algebra immediately 
yields, 



E = nu^^L^ + ^{l + g)N^ (24) 

Note that when L = and (7 = this is the standard result for TF energy of the ground 
state of a fermionic system in two dimension in an oscillator confinement. The effect of 
interaction is to scale the density and hence the energy by a factor which involves (1 + g) 
and the effect of rotations is contained in the L dependent term. 

The interaction energy alone may be separated, by computing the interaction term in 
the density functional using the TF solution for the spatial density. We have therefore, 

2N^g 

Eint = hujo — , (25) 
9^L^ + lN\l + g) 

Note that with the increase in angular momentum the effect of interaction decreases as it 
should. 

A few general remarks are in order here. In the actual many body system, there will be 
in general a tower of states in a given angular momentum sector. However, as is obvious 
from the expression for the total energy the TF method gives the lowest energy state or 
the ground state in each angular momentum sector. This is easily seen by taking the limit 
L >> N, when the energy is simply huJoL which is the energy obtained without any radial 
excitations. These states are often referred to as Yrast states in the literature. 

We briefly discuss the case with an arbitrary but very short range interaction. The leading 
term, as remarked before, is simply the contact term discussed above. The corrections 
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may be systematically worked out using the gradient expansion method |T5|,|T6|. Below, we 
mention the results from the next to leading order correction. Again the corrections may 
be obtained analytically. The results are given as corrections to the results obtained with 
contact interaction alone. 

The next higher order correction to the contact interaction in the gradient expansion of 



a short range potential is of the form ||T6 



VsR = -V2V^5{r) (26) 

where V2 is related to the second moment of the interaction potential. Note that in the 
gradient expansion only even moments appear in each term. The Thomas Fermi free energy 
functional is given by 

F[p] = J d\[^ + ImfirV(r) + ^p' - ^V^] -pj dhp{r) (27) 

Variation with respect to density gives the following self consistent equation for the 
density 

^2^ 1 
- V2V^p{r) + — - (1 + ^)p(f) + -m^]V2 = /i, (28) 

It can be easilly seen that for V2 = eq. ( ^8]) reduces to eq. (plQl) . 

The circularly symmetric solutions of the above inhomogeneous partial differential equa- 
tion can be found easily. Replacing V2 with the dimensionless coupling g' , where g' = ^^^3*^° 
we can write the solution as 

P = Pc + Per (29) 



where pc is given by eq.(21) with the contact term alone and pcr is the extra term which can 



be written in terms of modified Bessel functions 



2 uj^ lojkx) 



where 



9' 



(31) 



Here we have used the dimensionless variable x = f- with oscillator length scale L = as 

to o u mujo 

the unit of length. Here after we use the subscript " c" to denote the contribution from the 
contact term alone. 

The corrections to various other relevant quantities may be easily found with the density. 
The total no. of particles N can be similarly given by the following expression. 

The total angular momentum L ican also be written as 

The energy is again given by 



^ - ^ + II + 1^(1 - (34) 



The interaction energy is given by 



1 / ^ 

E,^t = E--{2 + ^{l- ^))cuL - hujo I d^rp\ (35) 



We note that the last term which involves the second moment of the density can be expressed 
in terms the modified Bessel functions. 



B. Logarithmic Interaction 

We now consider the case where two fermions interact via a logarithmic interaction. For 
quantum dot systems this is close to being realistic interaction at short distances where as 
the interaction is varies inversely as the distances at long distances [f^ . The two dimensional 



TF atom with logarithmic interactions was solved earlier by Bhaduri et. al. [^] and was 
applied to quantum dots recently by Sinha et. al.. 0) in order to explain the shell effect in 
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quantum dots. The extension of these results for the rotating case is quite staright forward. 
We give the main resuhs below: 
The TF mean field is given by 

1 I ~* ~^ I 

V{r) = -mn'^r^ - e) j ^Vp(r^)Zn( '"^ ~ ' ), (36) 

where a is an arbitrary parameter with the dimension of length. In order that the interaction 
is repulsive it is essential that the parameter a is much larger than any other length scale in 
the problem. The largest length scale in the problem is the turning point and we set a = tq. 
While this may seem arbitrary, in the analysis of the ground state of many fermions this 
choice is justified by a numerical analysis of the data on addition spectrum 0. The effective 
strength of the interaction is denoted by . Unlike in the three-dimensional Coulomb case, 
this parameter is not dimensionless, but has the dimensions of energy. The only change 
from the analysis of the ground state in the absence of rotations is in the first term where 
the effective frequency is different from the confinement frequency. We therefore only give 
essential results below and details may be looked up in Ref. 0. 

The spatial density is given by the solution of the differential equation 

^VV(r) = ^elpir) - ImV? (37) 
2m ' 

For a circularly symmetric spatial distribution, the above equation is easily solved 0. 
The solutions are given by the modified Bessel function as. 



1 /o(x) 
It (1 - ^) ^ h{x^) 



P(^) = ir7^-:JT[1-T^]> (38) 



where we have again introduced the dimension-less variables. 



with b = has the dimension of length 

with xo is again the classical turning point in units of h. 

The total particle number can then be obtained by integrating this density 

11 







^0 



where B = 

Since is fixed, the turning point Xq can be determined from this equation. Similarly 
one can calculate the total angular momentum. After some staright forward algebra it can 
be written as 

- = ^-[N{xl + 4) - -/-^(l - -^)] (41) 



The energy is given as 



E = /.- + [l + -^(l-^)]a;L (42) 



Therefore at a given N and L the constants ro and u are determined by the above set of 
equations. This is in general done numerically unlike in the case of contact interaction in 
order to express the energy as a function of N and L. The interaction energy also can be 
written in a closed form and is given by the following expression. 



IV. DISCUSSION OF RESULTS 

We have outlined an analytical procedure based on the simple Thomas-Fermi approach 
to the many body problem of rotating fermions in two dimensions. We have derived the 
density, the angular momentum as well the energy for model two body interactions. For 
simplicity we have chosen the contact interaction and logarithmic interaction. While these 
may not be the exact form of the interaction in the two dimensional systems they may under 
different circumstances approximately reflect the true nature of the interactions. 

We further illustrate these analytical results with numerical computation of quantities 
of interest by taking iV = 50 which is reasonably large for TF results to be valid. In Fig. 
1 we have plotted the change in angular momentum as a function of rotational frequency 
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measured in units of the confinement frequency. The curves represent the behaviour of L 
with the two model interactions. While in the case of contact interaction the relationship 
is almost linear all the way upto 0.6, the non-linear dependence of L on a; builds up much 
earlier in the case of logarithmic interaction. From eq.(^) and eq.(p2Dit is easy to see that 

^V(i - 4) oc LViv^ 

in the case of contact interaction. Hence the nonlinearity sets in only when the ratio between 
rotational frequency and the confinement frequency is considerable. 

As seen in Figs. 2 and 3, the total and interaction energy are sensitive to the behaviour 
of L as a function of u. We have shown the ratio of total energy with the same in the absence 
of rotations in Fig. 2. As expected the energy increases slowly at first, but becomes linear 
at high angular momentum since uL term starts dominating. This is a precursor of the 
expected behaviour in the Yrast region of the rotational spectrum which we have discussed 
earlier. It will be instructive at this point to mention the behaviour of the rotating nuclei 
in the Yrast region. There at very high angular momentum sector though the excitation 
energy is quite high, the nucleus remains cold as most of the energy is spent in rotating the 
nucleus. In the Yrast region the interaction is of lesser importance as can be seen from Fig. 
3. 

A related question of interest that is relevant to a rotating nuclei as well as the two 
dimensional fermionic system is that under heavy internal stress generated by the centrifugal 
and Coriolis forces generated in the Yrast region ( [|1^]) the density profile may drastically 
change. This issue has been investigated by Gallego and Das Gupta ( [|1^]) and requires the 
consideration of a non radially symmetric density. This question as well as system quantum 
corrections will be addressed in a forthcoming paper [|1^ . 
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FIGURES 

FIG. 1. Plot of the angular momentum versus the ratio —.The dotted line is the result for 

contact interaction, while the solid line gives the result for logarithmic interaction. The number of 
particles is kept same for all the cases at 50. The angular momentum is in units of h 

FIG. 2. Plot of the total energy E (scaled by £^ at L = 0) versus angular momentum L. Here 
also the solid line is for logaritmic interaction while the dotted line is for contact interaction. All 
other parameters are kept same as in Fig.l 

FIG. 3. Plot of the Interaction energy (Eint) again scaled by the Eint at L = as a function of 
the angular momentum (L). Here also the solid line is for logaritmic interaction while the dotted 
line is for contact interaction. All other parameters are kept same as in Fig.l 
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